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Various phenomenological models of particle multiplicity distributions are discussed using a gen- 
eral form of the grand canonical partition function. These phenomenological models include a 
wide range of varied processes such as coherent emission or Poisson processes, chaotic emission re- 
sulting in a negative binomial distribution, combinations of coherent and chaotic processes called 
OA , signal/noise distributions, and models based on field emission from Lorentzian line shapes leading to 

Lorentz/Catalan distributions. These specific cases can be written as special cases of a more general 
distribution. Using this grand canonical approach moments and cumulants, combinants, hierarchical 
structure, void scaling relations, KNO scaling features, clan variables and branching laws associated 
• with stochastic or ancestral variables are discussed. It is shown that just looking at the mean and 

fiuctuation of data is not enough to distinguish these distributions or the underlying mechanism. 
A generalization of the Poisson transform of a distribution and the Poissonian decomposition of it 
' into a compound or sequential process is also given. 
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^ . I. INTRODUCTION 

o 

' Pion multiplicity distribution and their associated fluctuation and correlations have been studied at the CERN SPS 
in the past, are being investigated at RHIC presently, and in the future will be studied at the LHC. Because of 
■ the very large number of pions and other produced particles in very high energy collisions, event-by-event studies can 
I be carried out and are important tools for many reasons. One current reason for studying them resides in the hope 
. that anomalous fluctuations will remain from a transition of a quark-gluon phase to hadronic phase, thus offering one 
,_H ' of the signals of the formation of the QG phase. Several models predict large fluctuations such as the disoriented 
^ , chiral condensate and in density fluctuations from droplets arising in a first order phase transition [^. A well 
2 ' known procedure for studying correlations uses the Bose-Einstein symmetries associated with pions in a Hanbury 
' Brown-Twiss analysis Such an analysis gives information about the space time history of the collision through 
^ measurements of source parameters. If the density of pions becomes large, Bose-Einstein correlation may also lead 
k>( to a strongly emitting system which has been called a pion laser @]. The pion laser model has been recently solved 
analytically by T. Csorgo and J. Zimanyi jsj. The importance of Bose-Einstein correlations has also been illustrated 
d in the observation of a condensation of atoms in a harmonic oscillator or laser trap Previous interest in pionic 
distributions have centered around the possibility of intermittency behavior and fractal structure based on parallels 
with turbulent flow in fluids. A distribution widely used to discuss such features has been the negative binomial (NB) 
distribution with its associated clan structure and KNO scaling feature |T^. KNO scaling properties 

have been interpreted in terms of a phase transition associated with a Feynman- Wilson gas Jl^. Various other issues 
associated with pions include evidence for thermalization p^, critical point fluctuations |17,18|, fluctuations from 



a first order phase transition |19[|, charge particle ratios and question of chemical equilibriunipOH , the behavior of 
fluctuations in net charge in a QG plasma for transition |2iy22t| . 

For lower energy heavy ion collisions, multifragmentation of nuclei takes place. The fragment distribution can also 
be described statistically by considering all the possible partition of A nuclcons into smaller clusters [^,|^ . This study 
give s a tool for the description of nuclear multifragmentation distributions psf , nuclear liquid-gas phase transition 
|^6| , critical exponent, intermittency, and chaotic behavior |p7| , p8| of nuclear multifragmentation. The same model 
can describe pionic distribution . This possibility arises in our approach which is based on Feynman path integral 
methods where symmetrization of bosons or anti-symmetrization of fermions leads to a cycle class decomposition of 
the permutations associated with these symmetries. The correspondance comes from the identification of clusters of 
size k and the cycles of length fc in a permutation as discussed below. 

In next section a summary of the generalized statistical model will be given. Various models of particle multiplicity 
distribution will then be developed. Moreover, we derive a generalized model which can further be reduced to a 
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geometric, negative binomial, and Lorentz/Catalan model. These various models used in describing pion and particle 
multiplicity distributions then are examined with the generalized model in Sect. III. In Sect. IV, we further compare 
various statistical properties between different models within the generalized model. 



II. GENERALIZED PROBABILITY DISTRIBUTION 



Consider a system composed of N different types of species or objects which could be the fragments in a 
fragmentation or in a cycle class description. Any event of such a system can be associated with a vector 
n = {uk} = {ni,n2, ■ ■ ■ ,Tik, ■ ■ ■ ,npf) or l"i2"^3"^ • • • fc"*" • • • A^"" where the non-negative integer Uk is the num- 
ber of individuals of species k. For example can be the number of clusters of size k or the number of cycles of 
length fc in a given permutation of n particles. The later is important for Bose-Einstein and Fermi-Dirac statistics. A 
general block picture of n is shown in Fig.|l|a. Fig.|l|b shows how the various partition can be developed as an evolution 
from successively smaller systems. The multiplicity, i.e., the total number of individuals is then 

N 

M = (1) 

k=l 

The number of species A^ can be infinity in general. 



□ 




FIG. 1. Building partitions with blocks. 



Various probability distributions related with this system can be developed by assigning an appropriate weight Xk 
to each individual of type k species. A weight W{x, n) is then given to each event n and the type of weight that will 
be considered has the structure; 



TV 



fc=i 



(2) 



The rife! are Gibbs factors. The Xk will be given below and contains various physical quantities. Summing the weight 
W{x, n) over all the possible events of n, the grand canonical partition function Z{x) of the system is given as 
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Z{x)^Y.^ix,n)=J2]l 



n fc=l 



exp 



■ N 
.fc=l 



(3) 



The last equation holds due to the form of Eq.(|^) of the weight W{x,n), i.e., the factor x'^'' /ukl is the n^-th order 
expansion term of e^* . 

Introducing other quantities ak to each individual entity or group of type k, a = {ak} — (ai, a2, • • • , ak, • • • , un), 
we can define a canonical partition function Za{x) with a fixed A as 
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^ = ^ akUk 

k=l 
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k=l 



(4) 
(5) 



with Z(x) — Za{x). Here X^a^ summation over all events with a fixed value of A, i.e., over a canonical 

ensemble of a fixed A, and the J^a ^ summation over all the possible values of A; it becomes X]a=o ^^^'^ '^^^'^ °^ 
ak = k with positive integer k. The case ak — k is encountered in fragmentation problems and permutation problems. 
If we take Xk oc 2"*=, then the canonical partition function Za{x) is the dependent term of the grand canonical 
partition function Z{x). The physics of the canonical ensemble depends on the choice of the quantity ak p3| , ^ . 
There always is at least one event having ^ = 0, i.e., the event where all n^'s are zero, n = 0. Thus Zo(x) = 1 if all 
ckfc are non zero positive since then there are no other possible events having A ^ 0. Due to the form of the weight 
W{x,n) given by Eq.(|2|) the canonical partition function Za{x) satisfies a recurrence relation [ p7|p9| 

Za{x) ^ -^^^akXkZA-akix) (6) 
k 

This relation is nothing but the constraint Eq.lQ) in terms of the mean < Uk >a using Eq.(^). For non-zero positive 
ttfc, there is no case having A <Q, i.e., Za — for A < Q and thus the Za can be obtained by the recurence relation 
of Eq. (^) starting from Zq {x) — 1. 

To study the statistical properties of a system where the number of individuals is independent of their species 
type, we can choose afc = 1. For fragmentation, all the clusters are treated the same independent of their size k 
or internal structure by choosing ak = 1. In pion distribution, 7r+, tt", and tt" are treated the same if we choose 
ak — 1 independent of their charge. By choosing afc = 1, we can study jet distributions without considering any 
further process of hadronization. Then A — '^^.Uk — M and the canonical ensemble is a set of systems having 
the same number of individuals A — M independent of their internal structure and the canonical partition function 
Za{x) = Zm{x) is the expansion of grand canonical partition function Z{x) in terms of z with a power of multiplicity 
M and z'^ — z*^ counts the total multiplicity of individuals in the system. As we will see later, this uniform treatment 
of all the species independent to their internal structure leads to a Poisson distribution. 

If we take au = k with the fc to be a positive integer representing the number of constituent particles (such as 
nucleons in a nuclear fragment or cycle length of a permutation), then A — X]fc ^"■fc the total number of the 
constituent particles in the system and z^ term in Z[x) is the canonical partition function Za with the total of A 
constituent particles. For this case the weight variable z is the variable counting the number of constituent particles, 
i.e., the fugacity of the constituent with z — and with /i the chemical potential of the constituent particle. This case 
is related to the nuclear multifragmentation model in Refs. p3| , p4| and various models of pion distribution discussed 
below. For ak — k, and if we also take Xk — x then all the clusters are treated the same independent of the size 
ak — k. The weight W oc a;*^ is then the same for all partitions having the same multiplicity except for the Gibb's 
factor. Thus x counts the multiplicity with x being a fugacity of a cluster. If we choose Xk = z^ then the weight 
W{x, n) is the same for all partitions having same A. All the nucleons are treated the same independent of the cluster 
it belongs to thus counting the number of constituent nucleons with z being a fugacity of the constituent. If we 
choose Xk = z"'' for a general value of ak then the weight W{x, n) is the same for all partitions having same A. All 
the constituents per unit ak are treated the same independent of the cluster it belongs to thus counting A as the 
corresponding total quantity of ak with z being a fugacity of unit ak- 

On the other hand, if we choose ak to be the energy of species k or fc-th level, then A = akrik = X^fc ^fc^fc = E 
is the total energy of the system and Za{x) becomes the partition function of a canonical ensemble with a fixed total 
energy E. If we choose ak to be the charge qk of species k, then A = J^k Qk^k — Q is the total charge of the system. 
Most discussions in this section applies to a general value of ak- However we will concentrate more on the choice of 
ak = k here for simplicity in notation. 



A. Probability distribution 

In a canonical ensemble of fixed A, we can define a probability distribution of a specific partition n as 

PA{x,n)= \' ' (7) 
Za{x) 
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With this probabihty, various mean values, fluctuations and correlations of the number of species Uk can be evaluated 
as a ratio of canonical partition functions for two different values of A such as [|3|j2^ 



rik'- 



(ufc - m)! (rij - ;)! 



E 



{uk - my. {uj - l)\ 



a, (^) 



Za{x) 



(8) 



Thus we have < n/c >a— Xk 



Za- 



j-g-^ — , which shows that the mean number of species /c in a canonical ensemble with 
fixed A is proportional to the weight Xk assigned to the species and the ratio of canonical partition function with 
different A. The recurrence relation of Eq.(||) then follows simply using the fact that A = J2k=i < nk >■ This 
distribution has been used in describing various fragment distributions in nuclear multifragmentation [psj-psf. 

Now knowing the partition functions Z{x) and Za{x), we can associate a probability Pa{x) of the system to have 
a fixed value of A in a grand canonical ensemble as 



Pa{x) 



1 



1 



Za{x) 

Z{x) ~ Z{x) r(A+ 1) 



-) Z{x,z) 



z=0 



Z(x, z) — Za{x)z'^ — exp 



E 



XkZ 



(9) 
(10) 



The last step follows from the fact that the power term of Z{x, z) is Za{x) if we put Xk = z"*". Thus the generating 
function Z{x, z) of Pa can also be looked at as a grand canonical partition function with the weight Xk replaced to 
be Xkz"'', where the variable z counts A explicitly and Z{x, 1) = Z(x). If we consider Pa{x, z) — Pa{x)z^ then z has 
two roles; one as a weight which is assigned the same to each constituent and another as a generating parameter of 
the probabihty PaIx). For the case that Zo{x) = 1, the void probabihty Pq is the inverse of the grand canonical 
partition function, i.e., 



Po{x) 



Zojx) 
Z{x) 



^Z-\x) 



(11) 



Inversely the canonical partition function Za{x) is the probability Pa normalized or rescaled by the void probability 
Po, i.e., 



ZAix) 



1 



TiA+1) 



Z{x, z) 



z=0 



Pa{x) 
Poix) 



(12) 



Another type of generating function of Pa that is frequently used may be defined as 

^ oo 

G{x, u)^Y. - = ^ E - 



Z{x, z ~ 1 — u) 



Z{x) 



exp 



Y^xk[{i-ur^i] 



Pa{x) 



T{A + l) 



(13) 
(14) 



We see that G{x, 0) = 1, G(f, 1) = Po(x) = Zq{x)/Z{x) and G{x, I - z) = Z{x, z)/Z{x). Once the probability Pa{x) 
is determined, various statistical quantities can be evaluated. Also the grand canonical partition function Z{x^ z) is 
given once we know the thermodynamic grand potential 



r2(af, z) = — In Z{x^ z) = — 



XkZ 



(15) 



This result can be used to study the statistical properties of a system. Moreover various moments and cumulants, 
mean values and fluctuations may be obtained using the generating function [E3^,e3. 
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B. Moments and cumulants; combinants and hierachical structure 



This subsection gives general expression for various quantities that will be used later when we discuss specific 
models. Since the probability of a specific event n in a grand canonical ensemble is given by 



P{x, n) 



W{x, n) 
Zix) 



(16) 



the mean of a quantity F in a grand canonical ensemble is related to the mean of F in a canonical ensemble as 



Za{x) 1 



^FW{2 



^ Pa{x) FPa{x, n) = J2 P^i^) <F>A 



(17) 



Thus once we obtain Za{x) using the recurrence relation of Eq.(^, then using the mean in canonical ensemble given 
by Eq.(H), we can obtain mean in grand canonical ensemble. Since, in a grand canonical enemble, A can be infinity, 
i.e., < ^ < oo with Za{x) — for yl < if all a/j > and —oo < A < oo if ak can be negative, we have 



Z{x) = ^ Za{x) = ^ ZA-mak-la, {x) 



(18) 



A A 

for integer m and I. Using this fact, we can see easily that 



<nk> ^nfeP(x, n) = ^PA(f) < Uk >a= ^-Pa(^) 



Xk- 



Za{x) 



<M> = 



Z{x) ~ 

n \ k / k 



Xk 



(19) 
(20) 



This result shows that the weight factor Xk in this model is the mean number < > in a grand canonical ensemble. 
The m-th power moment of A and its factorial moments are given simply by 



N 



<A'">(x) = ^^afcnJ P(a 



1 



r(^ + i) 

T{A-m+l) 



TjA+l) 
z^T{A-m + l 



-P{x, n) = 



Z(x) 
1 



Z-] Z{x:,z) 



z=l 



Z{x) Wdz 



— Z{x, z) 



z = l 



(21) 
(22) 



Similarly the m-th cumulants , which is the power moments of a^, and the factoral cumulants are 

/ ^ \ 



< a™ > (f ) 



/ N 



d 



dz 



In Z(x, z) 



/™(^)^ E 



r(Qfc + i) 

^ r{ak ~ m+1 



— — lnG'(a;, u) 
du ' 



We can see easily that, for the power moments of a^, 

N 



k=l 



Uk Xk 

2=1 k=l 



-rik ) = 



— In Z(x^ z) 
dz ' 



N 

E 



r(afe + 1) 



u=o r(afc -m+1) 



Xk 



(23) 



(24) 



(25) 
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< a > = ^ afcXfc = /i =< A > 
fe=i 

N 

<a^ > = ^ alxk = /2 + /i =< {A- < A > f >=< A^ > - < A ct^ 
fc=i 

N 

<a^> - ^a|xfc = /3 + 3/2 + /i (A- < A >)3 > 



(26) 
(27) 
(28) 



fc=i 



The power moments of ak are directly related to the power moments of A measured from the mean < A >, i.e., the 
cumulants < a"* > are same with the central moments of A. This simple relation does not hold for to > 4 but we 
can evaluate them starting from < a'^ > using the recurrence relation 



< > (f , z) 



dz 



m+l 



In Z{x, z) = [z-ij-] < a"" > {x, z) 



(29) 



Similarly the TO-th factorial cumulants /m, which is the factorial moments of a^, can be found using recurrence relation 



fm+ijx, z) 

fmi,'^-! z) 



^ ( ^ ) ln/„i(f, z)-m 



starting from 



N 



fo{x, z) = < M > (f , z) = In Z{x, z) = -Q.{x, -z) = ^ XkZ°'*' 



fc=i 



AT 



fi{x,z) = < A> {x,z) = Y^ 



akXkZ 



(30) 

(31) 
(32) 



fe=i 



The reduced factorial cumulants Km defined in Ref. pO| corresponds to the factorial cumulants /,„ normalized with 
mean number < ^ >= A as 



{x ^ z) 



/m(^7 z) 



(33) 



Thus with Ki — 1 and kq = fo =< M >. 

A Taylor expansion of {p + g)" with a real number a, w.r.t. p for q 7^ 0, is 



5;! 



^ ^ r(a+l) a-n n 

P=0 n=0 ^ ' 



E 

n=0 



(34) 



Thus, from the factorial structure of f„i and from Eqs.(|l3|) and (p4D, 



Z(a;, z = 1 - u) 

G(a;, -u = — — = exp 

Z{x) 



exp 



E 



{-uAy 



ml 



(35) 



Since Po{x) — G{x,l) the probability Pa can be obtained by 

A 



Pa{x) 



1 



r(A + i) 



-) Gix.u) 



T{A + 1) 



(36) 



when Km is independent of the mean < A >— A as used in Ref. Using the power moment < a™ > of Eq.(p3| 
we also have 



Zix) 



E 



< a™ > (x) 



(37) 
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Since /q = In Z =< a" >, 

Z{x, z) = G{x, u = l- z)Z{x) = G(f , u = l- z)e^°(^) 



exp 



E 



-/m(:E) 



exp 



.m=0 



(38) 



The generating function Z{x, z) differs from the generating function G{x, u) only by an extra term of m = in their 
exponent. 

The relation between the XkS and the Z shows that the XkS are also the combinants of Ref. |3^]. In turn the 
combinants can be related to the factorial cumulants defined by 



In Z{x^ z) = 



XkZ 



/ J I J" 



k m—0 

The factorial cumulants fm are the m-th order factorial moments of ak of Eq. 



|. Thus 



k=m ^ ' 



Xk 



(39) 



(40) 



for ak — k. The normalized factorial cumulants, i.e., the reduced cumurant, is 

= /™/^'" = (m - 1)!k^-i (41) 

for a negative binomial (NB) distribution. This result of Eq.(p|) shows that Km for NB has an hierarchical structure 
of a distribution at the reduced cumulant level which was realized for the NB distribution in Ref. |^^. This result 
will be generalized later. 

Also using the above power moments and factorial moments we can study voids and void scaling relation, hierarchical 
structure, combinant and cummulant properties which will be discussed below. 



C. Multi- fragmentation versus multiparticle production 

Since our approach was first used to discuss multifragmcntation and then later extended to include multiparticle 
production, we briefiy mention some of difference between multifragmentation and multiparticle production. 

In nuclear multifragmentation, afc = fc is the number of nucleons in a fragment and Uk is the number of fragments of 
size fc. The mean values of the total number of fragments Af = Uk and the total number of nucleons A = kuk 
are 

<M> = ^ < nfe ^ .Tfc = /o = InZ (42) 

k k 

< A> ^^k <nk>=^kxk^ h (43) 



In nuclear multifragmentation, the total number of nucleons A is usually fixed. In multiparticle production the A is the 
total number of produced particles and is not fixed. Thus in multifragmentation we consider canonical ensemble instead 
of a grand canonical ensemble. The mean multiplicity < M >a= 'Ylik '"^^ >A and < Uk >a— XkZ A-k{x) / Z a{x) in 
a canonical ensemble are related with grand canonical ensemble as 

oo 

< M > = ^<M>aPa (44) 
a=o 

CO 

<nk> ^ ^ < Uk >A Pa = Xk (45) 

A=Q 

The Xk =< nfc > in a grand canonical ensemble. But the weight factor Xk is different from the mean multiplicity of 
cluster size fc in a canonical ensemble < Uk >a- From Eq.(H), < Uk >a is the weight Xk multiplied by the ratio of 
ZA~k and Za, i-e., < Uk >a= XkZA-kjZA- In multifragmentation studies [p3|-p5[ we have used Xk — xz^ jk with the 
same weight x for each cluster and the same weight z for each nucleon. This choice of Xk gives Za{x, z) — ^^^-^^^ 
and thus ZA-k/ZA ~* as A — > oo where the z = e^^ with ^ the chemical potential. To determine the weight 
Xk experimentally, we should use a grand canonical ensemble, i.e., we need to consider various system with different 
values of A. 
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D. Clan parameters and void parameters and void scaling relations 



Van Hove and Giovannini have introduced clan variables Nc and ric to describe a general class of probability 
distributions, with most discussions of these variables centering around the negative binomial distribution |13| . These 
variables are defined as 

TV, = < M >= In Z = /o, n,=<A> /N,^ /1//0 (46) 

where the mean number of clans is Nc and the is the mean number of members per clan. The Z is the grand canonical 
generating function and thus Nc = J2k where Xk is the cycle class weight distribution x. The < A >= J2k k < rik > 
is the mean number of total members (particles). The Uk here is the number of clans of size k having k members. 

The clan variable Nc is also related to the void probability Pq = Zq/Z = 1/Z = e'^"; thus Nc = - \tlPq = fo. 
An important function in void analysis is x = — InPo/ < A >— Nc/ < A >— l/uc- Thus the void parameters, 
void probability Pq and void function v [ pO[ , are equivalent to the generalized clan parameters Nc and Uc with the 
equivalence given by: 

/o (x) = In Z{x) = - In Pq (x) = Nc (47) 
X(f)^i^ = ^-n,-i (48) 

The —X is the normalized grand potential — — /o for the mean < A>. 

Void analysis looks for scaling properties associated with x; specifically, x is a function of the combination < A > £^ 
where ^ is the coefficient of < ^ >^ in the fluctuation cr^ =< A > +^ < A >^. Since /2 =< > — < a >~< {A— < 
A > — < A > , the variance of ^ in a grand canonical ensemble becomes 

a"^ = < > - < A >^==< ^ alxk 

k 

= < A > +f2 =< A > +^ < A >2=< A>[1 + i{x) <A>] (49) 

with S, — K2, i-e., the normalized factorial cumulant. Since the variance for a Poissonian distribution is the same as 
the mean, cr^ =< A >, the ^ = 0; thus, the parameter ^ < ^ > represents a degree of departure from Poissonian 
fluctuation normalized by mean < A > of the distribution. A well known non-Poissonian example is a NB distribution 
which has C = ^ ^^'^ this becomes Plank distribution with x — 1. Using the recurrence relation Eq.(|30|) for /m, we 
can show that 

e(f , z)<A>ix)= x-\x, z) - z (^) lnx(f , z) - 1 = ^ ~ ^ ^ff - 1 (50) 

.,{x, z) = = ^A^{-d-z) < ^ > ] - 2^ (51) 

A NB distribution has x — lii(l + ^ < ^ >)/(^ < A >) while the Lorentz/Catala n (LC ) distribution discussed in Ref. 
Q and below has x = (^2^ < A> +1 - l)/(^ < A >). We will study in Sect. |IV A| , x ^ < A >, i.e., the void or 
clan variable vs the fluctuation for various choices of Xk summarized in Table |. 



E. Ancestral or evolutionary variables 



The LC model was shown to be a useful model for discussing an underlying splitting dynamics when ancestral 
or evolutionary variables p and fi were introduced into x and z as discussed in Ref. p3]. Percolation or splitting 
dynamics with a branching probability p and survival probability (1 — p) has a hierarchical topology as shown in 
Fig.^. Weighting each diagram by Xk = pCkp'^^^il — p)^ ^ the evolutionary or ancestral variables are related to the 
clan variables Nc =< M > and Uc =< A > /Nc- By taking Ck to be the number of diagrams of size k shown in FigH, 
the evolutionary dynamics is just that of the LC model. Then with j3 set equal to 1, xi = {1 — p), X2 ~ p{l — p) , 
X3 — 2p^(l ~p)^, X4 — 5p^{l — p)^, etc. The interpretation of this set of XkS reads as follows: xi has 1 surviving hne 
without a branch (p°(l — p)^) and one diagram (Ci = 1), X2 has 1 branch (p^) leading to 2 surviving lines ((1 - p)2) 
and one diagram (C2 = 1), X3 has 2 branch points (p^) leading 3 surviving hues ((1 — p)^) and two diagrams (C3 — 2), 
X4 has 3 branch points(p^), 4 surviving lines ((1 — p)^) and 5 diagrams (C4 = 5), etc. In these evolutionary /ancestral 
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variables the /o which determines Z is fo = J^^k = P for all p < 1/2. For p > 1/2, fo = J2^k is no longer a constant 
and is fo = Y^Xk = /3(1 -p)/p- 

Since the clan variables are Nc =< M >= fo and Uc =< A > /Nc = fi/ fo, then, for the LC model with evolutionary 
or ancestral varables, 



N, 



<M>=f3 



l-|l-2p| 



2p 



<A> =f3- 
1 



<A>/N,, 



2p(l-p) 



|l-2p|(l-|l-2p|) 



l-2p| 

It 



1 



2nc - 1 
(l±l)/2 



1 



(52) 
(53) 
(54) 
(55) 
(56) 



These can be reduced to Nc ^ P, Uc ^ (1 -p)/(l - 2p), and 2p = 1 - l/(2nc - 1) for p < 1/2 while Nc = P{l-p)/p, 
Uc = p/{2p — 1), and 2p = 1 + l/(2nc — 1) for p > 1/2. Since the branching probability p varies in < p < 1, the clan 
variable Uc has Uc > 1 wi th ric = 1 at p = and 1 and Uc = oo at p = 1/2. The behavior of Xk above p = 1/2 will 
be discussed in Sec. [II D| where a xq = 4>oo will be introduced. The LC model thus connects the clan variable Uc to 
the probability p of branching in the evolutionary or ancestral picture of Fig.||or in a percolation model. For Poisson 
processes p — Q (no branching ), Xk — f3Ski (only unit cycles and no BE correlations) and Uc = 1 (one member in 
each clan in average). 



1-p 
unit cycles 



period 2 cycles 



P 



P 



2p2(l -p)3 

period 3 cycles 




P 




5p^{i-py 

period 4 cycles 



FIG. 2. Evolutionary lines of descent in a hierarchical topology. Each branch increases the cycle length with probability p, 
survival 1—p- The probability distribution evolves from Poisson to chaotic. For clusters each branch generates a bigger cluster. 



F. Various distributions with Gauss hypergeometric series 

In describing either nuclear multifragmentation or multiparticle distributions, is taken to be a positive integer k. 
In a nuclear fragment distribution, an initial number of A nucleons are clustered into small nuclei of k nucleons. For 
multiparticle distributions, k can be related to the length of the cycle in a cycle class representation of permutations 
for Bosonic or Fermionic symmetry. Then A is the total number of produced particles of a given type such as pions 
and rifc is the number of cycle classes of cycle length k for a permutation of A identical particles. These cases have no 
situation with k — 0. If we consider a formation of M jets which are followed by a pion creation process, then k may 
be related to the number of pions from a jet. For this case, k — can also be included as a jet having no pion. 

Once we identify an appropriate Xk for a physical system, then we may use our general model to study the statistical 
behavior of the system. The various models used in pion distribution can be related to our general model with Uk ~ k 
by choosing Xk as a term in a Gauss hypergeometric series F{a, b; c; z); 



m=0 

r(a + m) r(a + n) F(a + rt + m — ?i) 



[a]n\a + n] 

m—n 

(58) 



r(a) r(a) r(a + n) 

Usefull values of [a]„ are [l/2]„ — (2n)!/(7'i!2^"), [1]„ ~ n\, [2]„ = {n + 1)!. Considering only positive k, we choose 
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T{a + k-l) r(c) r(6 + fc-l) 



[c]k-i (fc-1)! r(a) r(c + fc-l) r(5) (fc-l)! 

For this case, the thermodynamic grand potential or the generating function is 

oo 

fo{x) = fo{x, z) = log z) = -Q.{x, z) = ^ Xfc = xzF{a, b; c; z) 

If we allow jets without a pion, then we may allow k — also. For such a case, 

[a]k[b]k _ r(a + fc) r(c) T{h + k) z^ 
""^-'"^ [c]fe fc!"^ T{a) T{c + k) T{b) k\ 

OO 

/o(^) = /o(a:, 2) = log Z{x, z) = -ri(a;, z) = ^ Xfc = xi^(a, 6; c; z) 



(59) 
(60) 

(61) 
(62) 



We can see the only difference of the generating functions between above two cases is the extra factor z for the grand 
potential. We will mostly concentrate on the first case, i.e., fc 7^ 0. 
Using Eq. (|2^ ) or the recurrence relation Eq.(|30|) and Eqs.(|57|) - 

/„(a;,z) = -z™(^^^ n{x,z) = z^(^^^ log Z(x, z) = z™ (^^ ) xzFia,b; c; z) 

_ ^ \o\'m\b\m m+1 



z F{a + m, b + m; c + m; z) 



+xm z F[a + m-l,6 + m-l;c + m-l;z) 

The second order normalized factorial cumulant ^ = K2 and the void variable x then 

72(2;, 2) f2{x,z) 



(63) 



^{x, z) < A> 



<A> h{x,z) 

\Bl^z^F{a + 2, 6 + 2; c + 2; z) + 2^zi^(a + 1, 6 + 1; c + 1; z) 

^zF{a + 1, 6 + 1; c + 1; z) + F{a, b; c; z) 
Io{x, z) _ /o(.T,z) 
<A> fi{x,z) 

F{a, b; c; z) 



^zF{a + 1, 6 + 1; c + 1; z) + F(a, 6; c; z) 



(64) 



(65) 



For some values of a, b, and c, the hypergeometric function become a simple function; 

F{a,b;b]z) = (1 - z)"'' 
F(a,6;a;z) = (1 - z)"'' 

1 - (1 - zy-" 



F{a,l-2-z)^ 

z(l — a) 

F(l,l;2;z) = lim i^(a, 1; 2; z) 



ln(l - z) 



(66) 



Various models of pion distributions can be related with these functions as listed in Table ||. 



TABLE I. Various models with specific choice of Ofc = A: and Xk in hypergeometric series F{a, b- c; z) of Eq.(|^. Here fc = 
is not included, and thus fo —InZ = X^^^ Xk = xzF{a, b; c; z). Here 1 < < A'' with N —> ca except for Poisson which has a 
finite Nx. 



Model 


Xk 


Mx)=\nZ 


a 


b 


c 


Poisson (P) 


Nx5k,i or a; for fc = 1, 2, ■ • ■ , iV 


Nx = A 








Geometric (Geo) 


xz*" 


xz 
1-z 


a 


1 


a 
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Negative Binomial (NB) 
Signal/Noise (SN) 

Lorentz/Catalan (LC) 

Hypergeometric (HGa) 

Random Walk-Id (RWID) 

Random Walk-2d (RW2D) 
Generalized RWID (GRWID) 
Generalized RW2D (GRW2D) 



k, 



iy- 

-2(k-l) ( 2(fc-l) 

fc- 1 

fc! 

-2(k-l) ( 2(fc~ 1) 



2_2(fc-i) / 2(fc -^1) 

(fc-l)! "^^ 



a;2 

2 



(fc~l)! 



—a; ln(l — z) 


1 


1 


2 


2x[l - (1 - z)^^^] 


2 


1 


2 




d 


1 


2 




1 
2 




1 (h\ 


xzF{\,\;l-z) 


1 
2 


1 

2 


1 


xzil-z)'" 


a 


b 


6 


xzF(a, a; 1; z) 


a 


a 


1 



TABLE II. Factorial cumulants for various choices of of Table^ 



/o = log Z 



fi=<A> 



Nx = A 

— a;ln(l — z) 
jK^- _xln(l -2) 
2a;[l - (1 - 

7- 2 

■^(l-z)l/2 



a; 



x-^ 

1 — 2 

T 2 

■^(l-z)l/2 



a; 



(l-z)» 

+ -- 



2 (l-z)l/2 2 (1-2)3/2 

T (a-l)z I ^ az 
(l-z)° 



'(l-z)3 



^ (1-Z)2 

J' (1-Z)3 + ^(l_^z)2 

•^(l-z)3/2 
^ a 



4 (l_z)3/2 4-^(l_z)V2 
^ (a-l)a I _£i(o + l)_ 



for m > 2 

Xn-z)- 

(l_z)-. l^ci-z) +^ 

•^(l-z)'"~l/2 
_[-l/21™ ^ [l/2]„ 
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More detail discussions and related physical systems of these distributions will be given and discussed in the next 
section. Brief discussion about the weight Xk of each model follows. A Poisson (P) distribution is generated when 
there is monomers only, i.e., Xk = Nx5k,i which gives /o = Nx. A Poisson can also be generated if all the clusters or 
cycles are treated the same with the same weight oi xt — x independent of their size fc for 1 < fc < which also gives 
/o = Nx. The second case can also be viewed as weighting each constituents by x'^/^ so that Xk — {x^^^)^ — x. The 
weight X counts the multiplicity M and the multiplicity of clusters (or of monomers for Xk = NxSk,i) has the Poisson 
distribution. All the other distributions considered in Table | have several factors in the weight Xk- One factor in 
Xk is z'^ which is a fc dependent geometric term and comes from assigning the same weight z to each constituents 
independent of the cluster or cycle classes it belongs to. Another factor of weight is independent of k such as the x 
in Table | which comes from assigning the same weight x to each cluster or the cycle class as a whole independent of 
its internal structure. These two factors, xz'', are multiplied by a fc dependent or independent prefactor. A geometric 
(Geo) distribution follows when there is no other weight factor beside x and z, i.e., no k dependent prefactor so that 
Xk = xz^ . The Geo with z = 1 for a finite N is the same as the Poisson distribution; both have /o = Nx. The 
negative binomial (NB) which appears frequently in various studies has a weight factor assigned to a cluster or cycle 
class given by xz^ /k. This has an extra size dependent factor of 1/fc compared to the geometric distribution. The 
signal/noise model (S N) ha s a two part structure and interpolates between a Poisson and NB distribution as will 
be discussed in Sect. [II C| . The geometric distribution is the signal component of SN while the NB distribution is 
the noise component of SN. The Lorentz/Catalan model (LC) has in its weight a shifted Catalan number divided by 

22(/c-i)^ that is — ^ '^^^ ^' ^ ^ , beside the xz'^ factor which is the weight for Geo model. The Catalan 

numbers given by i 1 I {k + 1) are 1, 2, 5, 14, • • • for fc = 1, 2, 3, 4, • • • and the shifted Catalan numbers given by 



^/f 1^ ) ^'"^ ^' ^'^^^ ' ' ' '^'^^ importance of this factor is shown in Fig. |^ of section II E. As can be seen 

from the arguments of the hypergeometric function F{a, b; c; z) in Table ^ the hypergeometric model with b — 1 and 
c = 2 (HGa) include Geo, NB, SN, LC as a special case of HGa depending on the value of a. Other models listed 
in Table || are based on random walks. The use of random walk results was originally due to Feynman [^3| in his 
description of the phase transition in liquid helium. The random walk aspects arise when considering the closing of 
cycle of length fc. We include them for completeness. Since the random walk in 1-dimension (RWID) is the same as 
LC except the missing 1/fc dependence compared to LC, RWID can be extended to a generalized RWID (GRWID) 
similar to the generalization of LC to HGa. A random walk model in 2-dimension has an extra factor of a shifted 
Catalan number and fc2~^*^''~^) factor compared to RWID and can also be generalized to GRW2D. 

Since fc = is excluded here, the partition function for these models are given simply by a hypergeometric function 
as Z = exp[^^.^ Xk] = ex.p[xzF{a, 6; c; z)] with various choices of a, b, c. For example the LC model has /o ~ 
Sfc^fc = xzF(l/2, 1; 2; z) and the NB model has /o = xzF{l,l;2; z). The geometric model Xk = yz'^ has /o — 
yzF{a, 1; a; z) = yzF{2, 1; 2; z) while the SN model is a combination of the geometric plus NB cases. These functions 
are special cases of /o = xzF{a, 1; 2; z) of the HGa model. The generalized random work in 1-dimension (GRWID) 
has /o = xzF{a,b;b; z) and the generaized RW in 2-dimcnsion (GRW2D) has /o — xzF{a,a;l; z). The factorial 
cumulants fm for these models are summarized in Table ||. The cases with c = 2 have a canonical partition function 
Zn which can be writen in terms of confluent hypergeometric functions U{u,v;w) and standard factor z"/ri! |34|]. 



G. Generalized model of Hypergeometric (HGa) 



We consider in more detail the hypergeometric model with b = 1 and c — 2 (HGa) here since it includes the NB, 
Geo, and LC models as special cases. This generalized model is related with the hypergeometric function with b = 1 
and c = 2 with an arbitrary value of a, i.e., F(a, 1; 2; z), and has the weight of 



_^^fc [fl]fc-i ^^z^T(a + fc) 



fc! 



fc! r(a) 



Its associated grand canonical partition function 



Z{x, z) — e 

is shown in Table |. From Table ||, we have 



exp 



.(a-1) V(l-z)(«-i) 



(67) 



(68) 
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r(a + m- 1) 



a+rn— 1 



r(a) (1 - 

^ ' ^ ^ ' r(a) V 2; y 



r(a + TO - 1) K2'"^(a;,z) 



r(a) 



(69) 



(70) 



The normalized factorial cumulant, i.e., the reduced cumulant shows the hierarchical structure of HGa at the 
reduced cumulant level with Am — a~*^™~^^r(a + m — l)/r(a) where Km is related to K2- This property was realized 
for the NB distribution in Rcf. |Q which is obtained for Eq.([70|) with a — 1 giving Am — {m — 1)1. The result of 
Eq.([70|) is a generalization of the NB result. 



Some moments for HGa are 

<A>{x,z) = fi{z,x) : 



xz 



Xix,z) 



(1 - z)" 

1 (l-z) 



C{X,Z) = K2{X,Z) 



fi{x,z) (1-a) z 
f2{x,z) a 



f1{x,z) X 



Since these relations give 



^(x, z) < A> [x, z) — K2(x, z) < A> {x, z) 



[{l-zf-^-l] 

(1 -z)("-l) 

f2{x,z) az 



(71) 
(72) 
(73) 

(74) 



fi{x,z) (l-z) 

the void parameters can be obtained in terms of the normalized fluctuation ^ and the mean number < A >— A hy 

_ f2/A _ ,f2 



<A0 = 



a + ^A a + f2/A /z + aA 

(a-l) 



z ^ \a + ^A 



aA^ 



aA 



h \aA + f2 



a-l 



/o(A,e) = iogZ(Ae) 
1 



X 



>im{A, ^) 



k 

A 



a-l 
a 



1 + — 
a 



[l-o.)iA 



1 + ^ 

a 



- 1 



a-l 

1 aA 
(T^X 



A 

aA 



1+4 

aA 



- 1 



fm{A,0 ^ T{a + m~l) n 
r(a) \a 



A" 



(75) 
(76) 
(77) 

(78) 
(79) 



for a given mean value of < ^ >= A and the fluctuation ^ or ^A or /z. 

Table | shows that the generalized HGa model becomes the Lorentz/Catalan (LC) model with a = 1/2, the negative 
binomial (NB) model with a — 1, and geometric (Geo) distribution with a = 2. However the NB should be considered 
as a a — > 1 limit of HGa; 



lim /o = -.Tlog(l - z) 
lim Z = (1 - zy 



limfo{A,0^x\n{l+^A)^x\n{l + f2/A) 
limx(Ae) = -^log(l + e^) 



(80) 
(81) 
(82) 

(83) 



Thus a NB distribution has x = ln(l + ^A)/(CA) with x = 1/^ while the LC distribution has x = (vl + 2^- 1)/(^^) 
with X = a/1 + 2^-4/(2^) and a Geo distribution has y = 2[1 - (1 + CA/2)-^]/{^A) with 1/x ^ ^(1 + ^A/2)/2 as 
limiting expressions of a more general x given by Eq.dTq) for HGa. 
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III. PION PROBABILITY DISTRIBUTION 



In this section we discuss various features of the pion multipilcity distribution and its associated fluctuations. In the 
introduction we noted the importance of such studies. The present work is a continuation of the approach developed 
in Ref. |^ which is based on Feynman path integral methods ||3^ . These methods lead to a cycle class decomposition 
of any permutation, with permutations appearing when Bose-Einstein or Fermi-Dirac symmetries are included into 
the density matrix. Any permutation of a particle has associated with it a vector n of which is the number of 
cycles of length k. A cycle of length 2 would be 1 ^ 4 ^ 1, etc. The constraint A = kuk must be satisfied for a 
canonical system. The M — Uk is the multiplicity which varies from 1 to A. A weight WAin, x) of Eq.(||) is then 
given to each n with the Xk in x the weight assigned to a cycle of length k. The grand canonical partition function 
Z{x) is Eq.(|^) and the probability distribution Pn{x) of n = J^k ^"■fc pions is Eq.@. The mean number < n > and 
the fluctuation of pions can be related to the Xfc's as 

< n > — ^ kxk 

k 

= =< > - <n >^=< n > +^ < n >^= ^ k^Xk (84) 

k 

from Eqs.@ and 



A. Poisson distribution and cycle of length 2 correlations 

A widely used distribution which also is a basis for comparison is the Poisson distribution which is obtained for 
the XkS by either having unit cycles only, i.e., Xk = xiSki or uniformly distributed cycles independent of their cycle 
length which is equivalent to the weig ht [xi/Ny/^ of each pion in a cycle class with size k (see Table ^. For this case 
n = M and the cycle class multiplicity M has a Poisson distribution. Then < n >= xi, =< n >, and P„ is 

Pn =< n >" e~<">/n! (85) 

with Z„ = Xi/n\ and Z — e^^. The probability of no events, the void probability distribution has n = and 
Pq = e~^"^. Poisson distributions appear in coherent state emission and are the Maxwell-Boltzmann limit of 
Bose-Einstein and Fermi-Dirac distribution in statistical physics. Fermi-Dirac and Bose-Einstein statistics lead to 
departures from Poisson results. The voids parameters are ^ = and x = 1 for this distribution. 

The first correction to Poisson comes from cycles of length 2. Then from Eq.(^4|) < n >= xi + 2x2 and — 
xi -t- 4x2 7^< n >. The P„ distribution can be generated for x = {xi,X2, 0, 0, • • •) using the recurrence relation Eq.(H) 
to obtain Z„ and Z(x) — exp[a:i -|- 2:2]. The Z„ is simply 

Zn = -(xiZ„_i + 2a;2.^„-2) (86) 
n 

The Z{x) — ^„ Zn becomes the Hermit polynomial e^^^~^' = J2n ^,^nix)z" when xi = 2xz and X2 = —z^. 

A system of one type only gives a Poisson distribution and a system of two or more types exhibits non-Poissonian 
behavior. Departures from Poisson statistics have been noted in level densities where random matrix theory gives the 
Wigner distribution which is not a Poisson distribution. The next subsection discusses a simple distribution which is 
the negative binomial and it can have large non-Poissonian fiuctuations. 



B. Negative binomial (NB) distribution 

A distribution frequently discussed for situation that depart from Poisson statistics is the negative binomial (NB) 
disribution. The cycle class vector x which generates the NB distribution is Xk = xz^ /k, logarithmic distribution, and 
Xk contains cycle lengths k of all sizes. For the choice (see Tables | and Z ~ [1 — z)^^, Z„ = ^T(x + n)/r{x), 
< n >= xz/il — z), — xz/{l — z)^ =< n > (1+ < n > /x). The P„ can be written in the well-known form 

z" T{x + n) _fn + x-l\f 1 Vf <n>/x \" 
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The X = 1 limit is also referred to as a Planck distribution since the fluctuation iscr^ =< n > (1+ < n >) and has the 
characteristic Planck behavior. Departures from Poisson statistics arise from chaotic sources as discussed in quantum 
optics. Intermittency and its association with chaos have also been studied using the negative binomial distribution 
pof . The voids parameters are ^ = j and x = '"^ "^^ ^'^^ ^^^^ distribution. 

Ref. gives several sources for the origin of the negative binomial distribution. These sources include sequential 
processes that arise from a composite of logarithmic and Poisson distributions, self-similar cascade processes and 
connections with Cantor sets and fractal structure, generalization of the Plank distribution, solutions to stochastic 
differential equations. Becattini etal have shown that the NB distribution arises from decaying resonances. The 
a-model of Ref. |lO), which is a self-similar random cascading process, leads to NB like behavior. The stochastics 
aspects of the NB distribution have been discussed by R. Hwa . Hegyi |Q has discussed the NB distribution in 
terms of combinants. 



C. Signal/noise (SN) model; Coherent and chaotic emission 

The signal/noise model (SN in Table | and ||) offers a continuous connection between a Poisson distribution (x = 
with y ^ and z ^ 1 keeping < n >— yz/{\ — z)^ finite) and the Planck or Bose-Einstein (BE) (x = 1 and y — Q) 
or more generally a negative binomial distribution (arbitrary x with y = Q). Since a coherent state is a Poisson 
emitter and a Planck or NB distribution comes from a chaotic state, the SN model interpolates between coherent 
and chaotic emission. The SN model is discussed in Ref. |^ and its connection to the cycle class picture developed 
here was initially shown in Ref. |Q. Here, further developments of it are discussed. To begin with the cycle weights 
have a geometric piece yz'^, distributed with z'^ , and a NB piece xz^ jk or Xfc = (y -I- {x/k))z'^. The y and z can be 
redefined in terms of a coherent signal (Poisson emitter) variable S and a thermal Bose-Einstein noise variable N 
using y = {n/x)(i+n/x) ^^"^ ^ ^ {i+n/x) ■ "^^^ probabihty distribution is 



-s 



l + N/x 



'^iN/x)il + N/x)^ ^^^^ 



with X = 1 being the Glauber-Lach model and L° is an associated Laguerre function. When — > the Poisson limit 
is realized with Xk = S5k,i and when — > 0, the distribution goes into a NB distribution. In terms of S and N the 
mean 

<n> ^ + ^S + N (89) 



while the fluctuation 



= < n > +^ < n >^ 



{l-zf 

yz xz xz^ 2yz^ 



(1-Z)2 1-Z (1-Z)2 (1-Z)3 



^ + iV+^ + ^ (90) 



X 



The void parameters are 



h ^ N {2S + N) 
fl ^ X [S + NY 



(91) 



^ io ^ x\n(l + N/x) + S/(l + N/x) 
^ h S+N ^ ' 

for this distribution. The SN model has important application to quantum optics and, in particluar, to photon counts 
from lasers Biyajima has suggested using it for particle multiplicity distribution as does Ref. [^. When the 
noise level ^ 0, ^ — > and x ^ 1- When the signal level S* ^ 0, f ^ l/x and X ^ ln(l -I- N/x). 

D. Lorentz/Catalan (LC) distribution; underlying splitting probability 

Probability distribution associated from field emission from line shapes which are Lorentzian have appeared in quan- 
tum optics. Its cycle class description was developed in Ref. Q and is Xk = xCkz'' /2'^^''''^^ with = ^ "^^^ _ j /k 
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a shifted Catalan number, distributed with NB with x = 1/2 (see Eq.(108) in Sect. [VC). The partition function is 
Z{x, z) = exp(2a;[l — (1 — z)^/^]) and the probabiUty distribution associated with it is 



Pn{x, z) = e-2-(i-^^)(2x)2"f/(n, 2n- Ax)^ (93) 

TV. 

with C/(m, v\ w) a confluent hypergeometric function. The void parameters are ^ — 2x^1- z ^^'^ ~ 
[VI + 2^ < n > - 1] with z = 2^ < n > /(I + 2^ < n >) for this distribution. 
The results of quantum optics, in the notation of Ref. can be obtained j32) when x = Tn/2, z = IW^jV? , 
= 7^ + 2W'-f. The W is an integral of the Lorentzian line shape r(w) = b/[{u! — wg)^ + 7^]i T is the time, and 
2xVT^ = iT. 

An interesting feature of the model arises when evolutionary or ancest ral v ariables are introduced which contain a 
branching probability p, survival probability (1 — p) as discussed in Sect. [IE, Percolation or splitting dynamics with 
evolutionary parameters p and /3 become equivalent to the LC model with clan variables Nc — (3 and Uc = li-2p) 

for p < 1/2. Specifically x — P/'lp and z — Ap{l — p) giving Xk — pCkp''^^{l — pY ■ In these evolutionary/ancestral 
variables the fo = = 2x{l — ^/l — z) which determines Z has the interesting property fo = ^Xk = P for all 

p < 1/2. For p > 1/2, — ^Xk is no longer a constant and is — ^Xk — /5(1 — p)/p. To keep /o = X^fe ^fe ^ 
constant without changing <rt>=Xj.fc<?ife>= [3 \i^2p\ ' ^ ~ introduced in Ref. The (poa = for 

p < 1/2 and (poo — f3{2p — l)/p for p > 1/2. For p > 1/2 there is a finite probability that the splitting will go on 
forever. In percolation above a certain p an infinite cluster is formed and (poo is similar to the strength of the infinite 
cluster. Moreover the sudden appearance of (poo is similar to the behavior of an order parameter in a phase transition. 
The appearance of (poo can also be interpreted as the sudden appearance of a jet without pion (fc = 0), when xq = (poo ■ 
A connection of the LC model can also be made with a Ginzberg-Landau theory of phase transitions and a Feynman- 
Wilson Gas. These connections are discussed in Ref. 13411. 



E. Thermal models 



Thermal emission of pions based on statistical mechanics and equilibrium ideas have been popular descriptions of pi- 
ous coming from relativistic heavy ion collisions. For thermal models the Xk = {VT^ /2TT'^){m/T)'^K2{km/T)/k for 
a cycle length fc or a cluster of size k with K2 a Mac Donald function. For low temperatures, Xk = {V/\'^)e~^"^/'^ /k^^"^ 
and the Boltzmann factor in mass, e^*"'™/^, suppresses large fiuctuations. In the high temperature limit and/or zero 
pion mass limit Xk — {V / 'k'^)T^ / k"^ . The Xk can be used to generate the pion probability distribution f „. The thermal 
models can be combined with hydrodynamic descriptions and an application was given ]3^ to 158A GeV Pb+Pb 
data measured by the CERN/NA44 and CERN/NA49 collaborations. The results of Ref. [^sj showed a Gaussian 
distribution with a width about 20 % larger than the Poisson result. For Table II of Ref. has xi = 260, X2 = 9.957, 
X3 = 1.097, 2:4 = 0.163, • • •. The cc^-'s have a rather sharp fall off with k and two terms xi and X2 give a reasonably 
good description of the Gaussian. A Poisson emitter with large < n > has also a Gaussian distribution with a width 
such that n >. For the first two x^'s just given xi = 260, X2 = 10, the J2^k — 270, < n >= ^ kxk = 280 and 

= 300. 



F. Fisher exponent r and other distrubutions 



Various forms for Xk encountered for the cycle class weights have a form Xk — xz^ /k'^ . The NB has r = 1, 
the geometric has r = and the SN has a combination. The LC model has r = 3/2 for large k using Stirling 
approximation. In clusters yields r is called Fisher critical exponent and it determines the droplet yields < Uk > 
around a critical point of a first order liquid-gas phase transition. It has been studied extensively in medium energy 
heavy ion collisions |42[ . Here we discuss its importance in particle multiplicity distributions. A power law of < > is 
the power law of Xk according to Eq. ( |l9|) in grand canonical ensemble and Eq. (|^) in canonical ensemble with finite size 
effect through Za- The r determines the grand canonical partition function Z{x) — cxp[Xj,a:fc] = exp[a;Xj, z'^/fc'^]. 
Bose-Einstein condensation of atoms in a box of sides L of dimensions d have x = L'^/X^ with \t = h/ {2'KmkBTy^'^ 
and have t = 1 + d/2. Feynman used random walk arguments, the closing of a cycle parallels a closed random walk, 
to discuss his choice of Xk in his discussion of a superfluid phase transition in liquid equilibrium. We have therefore 
consider two other case for Xk which asymptotically have r = 1/2 and t = 1. These are called random walk (RW) 
in ID and 2D respectively. The RWID has an Xk which is just k times the Xk of the LC model. The RW2D has the 
same r as the NB distribution in the asymptotic limit of large k. 
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G. Pion laser model 



A simple emitting source model for pions was introduced by Pratt |Q and solved analytically by Csorgo and Zimany 
1^. The model can be used to study the role of Bose- Einstein symmetrization effects on pionic distributions. An 
intersting property of this model is that a Poisson emitter of pions with strength ij can become a pion "laser". 
A manifestation of this property is a large enhancement of pions in a zero momentum mode. For this model the 
Xk = (?7/ryc)Vfc/(l - (7-/7+)''^^)^ with 7± = i(l + y ± V^TTy), y = 2R'^mT - 1/2 and r,c = 7+^^ For very large 
k, Xk = /k 1^^, a NB behavior with a; = 1 or the Planck limit. The z = rj/ijc- For small k and typical values of 
R and T the Xk = {R'^mTf/'^{'n/{R^viTfl'^Y/k'^. The l/k"^ dependence also appears in thermal models with zero 
mass pions. 



IV. COMPARISON OF PROBABILITY DISTRIBUTIONS WITHIN HGA 

The hypergeometric case (HGa) which we studied more in detail in Sect. II G included various cases of Table |; 
Geo, NB, LC which are distinguished by one parameter a of the hypergeometric model. Thus we use HGa to compare 
various models for pion distribution and other particle distributions. 



A. Voids and void scaling relations 

Void analysis looks for scaling propoerties associated with x = /□/ < n >; specifically, x is a function of the 
combination ^<n>=/2/ <n> where ^ = K2 is the coefficient of < n >^ in the fluctuation =< > — < 
n >^=< n > +^ < n >^. For example, a NB distribution has x = + £, < n >)/(f < n >) while the LG 
distribution has x = (\/2C < n > +1 — l)/(^ < n >). These two cases are limiting expressions of a more general x 
given by 

= (l + e<n>/ar°-l 
^ il-a)^<n>/a ^ ' 

of the HGa model. The NB x is obtained from the limit a — > 1 while the LC x follows for a = 1/2. 

We show the void variable x as a function of ^ < n > in Fig.^. This shows that we can vary a to fit data. Ref. 
pO| claims the NB distribution fits reasonably well the void distribution for single jet events in e'^e^ annihilation but 
Fig.^ shows all of the curves might fit such data up to ^ < n >~ 1 since the various curves are reasonably close up 
to this ^ < n >. Thus further investigations of this data is required to distinguish various models. Higher moments 
or cumulants might need to be compared for this purpose. Within HGa, K3(n,f,a) = (^^) for a given value of 
< n >= n and ^ according to Eq.@. Thus A3 = ks/C^ = (a + l)/a = 3, 2, 3/2, 4/3, 5/4 for a = 1/2, 1, 2, 3, 4 
independent of ^. The K3 may help in distinguishing various models for the data with a small value of ^, i.e., the 
region of < 1 in Fig.|^. The differences of Am between models with different value of a becomes larger as the order 
m of the reduced factorial cummulant becomes higher. 
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f<n> 

FIG. 3. X vs ^ < n > for a = 1/2 (solid line; LC), 1 (dash; NB), 2 (dash-dot; Geo), 3 (dot), 4 (dash-dot-dot-dot). For 
Poisson distribution ^ = and x — 1- 



B. Probability distribution 



Once we know the generating function Z{x) or fo{x) we may study the probabiUty distribution P„ using Eq.(|^) or 
the recurrence relation Eq.(a); 



P„ix,z) 



Z{x, z) Z{x, z) n! 



— ) Z{x,z) 



(95) 



for n — fcnfc. For NB which is a special case of HGa in the a 1 limit, Xk = xjr, 



7NB 



{x,z) = {l-z)-^ = P^\x,z) 



1 z"r(x + n) 
Z^ T{x) 



= ii-zy 



,z" r{x + n) 

^ r{x) 



(96) 
(97) 



The Pn for various cases of Table || are shown in Fig.|| and Fig|| with the same mean value < n >= n and the 
same fluctuation ^ = K2 = f2/n^ for fixed < n >= 10. Fig.|| shows KNO plots of < n > P„ versus n/ < n > for 
fixed < n >— 10 and 20. Fig.^ shows that the various models considered here have almost the same distribution 
for small fluctuation (^ = 0.01) and in this case they are very similar to a Poisson's distribution. For ^ = 0.05 the 
models are similar to each other except for larger n/n even though they are different from a Poisson distribution. For 
larger fluctuations such as ^ = 0.5, the models have very different forms even if they have the same mean value and 
fluctuation. Fig|| shows that the probability distribution of these models differ in their form for fluctuations larger 
than ^ « 0.2. 
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n/<n> 

FIG. 4. P„ for fixed < n >= 10 for a = 1/2 (solid line), 1 (dash), 2 (dash-dot), 3 (dot), 4 (dash-dot-dot-dot) and for ^ = 0.01, 
0.05, and 0.5 in log scale on the left and linear scale on the right. For ^ = 0.01 all distributions become very close to Poisson 
(thin solid curve) already. (Po becomes large for large ^.) 
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FIG. 6. KNO plot of < n > P,, for fixed < n 10 and 20 for a = 1/2, 1, 2, 3, 4. 

The KNO behavior for n — > oo can also be studied for the models listed in Table ||. According to KNO scaling, for 
distributions with a large mean value, the distribution < n > P„ becomes model independent in the new variable n/n, 
i.e., variable scaled by mean value. In general, any distribution becomes Gaussian for large mean < n > according to 
the central limit theorem; specifically 



Pn{n,a) 



1 



V2t 



: exp 



(n — n)^ 



1 



-1/2 



2a2 
1 



: exp 



1 



(98) 



with the mean < n >= h and the standard deviation a which is related to the reduced factorial cumulant f by 
—< > — < n n + ^n^. This means that the KNO scaling follows when the fluctuation is given by 

cr^ = £_n^ with a constant ^ or when ^ 1. Thus to compare KNO scaling properties of nPn for different models, 
the fluctuation of these models should have the same value of ^ = For a small ^, the Poisson component of the 

fluctuation cr^ = n + S^tt? becomes dominant and thus the KNO scaling behavior is broken. For large ^ the Poisson 
component is negligible and KNO scaling is realized, i.e.. 



1 



: exp 



1 /n \^ 



(99) 



The KNO plot of Fig.y shows that the various distributions have no KNO scaling property for small fluctuation 
(^ = 0.1) but show a KNO scaling property for large fluctuation (^ = 0.5). The effect of mean on nP„ is larger than 
the difference between different model for ^ = 0.1 while the effect of mean on nP„ is much smaller than the difference 
between different models exhibiting KNO scaling behavior for ^ = 0.5. Ref. shows the total charge distribution 
in hadronic collisions. From their fits we can extract the corresponding fluctuation which are ^ = 0.05 0.5 and 
n = 6 ~ 13. This means that the KNO scaling behavior is marginal for these data, i.e., just fitting the distribution 
nPn of the data does not show clear evidence for KNO scaling. We need to evaluate the explicit values of the mean 
number n and the fiuctuation ^ to check the KNO behavior of this data; the value of £,h should be large enough to 
show KNO scaling behavior as can be seen from Eqs.(|9^) and (p9|). 
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C. Sequential procedures and compound Poisson distributions 



A Poisson distribution plays a very important role in physics. As already noted, in statistical physics, Maxwell- 
Boltzmann statistics leads to Poisson probabilities. Other distributions are compared to the Poisson distribution 
which acts as a benchmark for comparison. The distributions considered in this paper can have large non-Poissonian 
fluctuations. The purpose of this section is to show how they can be rewritten as a compound process or sequential 
process involving one aspect that has a Poisson character. As an example the negative binomial distribution can be 
obtained from a compound Poisson- logarithmic distribution as discussed in Ref. Here, we extended this result to 
include the other distributions and we also show that the final distribution can be obtained from compounding it with 
another distribution, such as the negative binomial. In general, the underlying picture for a sequential process involves 
a two step procedure in which the observed particles arise from the production of "clusters" with the subsequent decay 
of each cluster producing its distribution of particles. The final distribution is obtained by compounding the probability 
distribution of the clusters with another distribution coming from each cluster and suming over clusters. 

The LC model was shown to be a useful model for discussing an underlying splitt ing dyna mics w hen ancestral or 



evolutionary variables p and jS where introduced into x and z as shown in Sects. II E and [II D and discussed in 
Ref. 1^ . The LC model thus connects the average size of clan Uc to the probability p of branching in the ancestral 
picture of Figj^. For p < 1/2, the mean number of clans is Nc = (3 and the mean number of members per clan is 
ric = {1 — p) / {I — 2p). For Poisson processes p = 0, Xk = (35ki (only unit cycles and no BE correlations) and ric = 1. 

Further discussions of the mean number of members per clan or of will now be given. The Xk =< > with 
< > the mean number of cycles or correlation of size k (c.f., Eq.(p^). The ti/j can also be considered as the number 
of clusters of size fc. This is a natural identification when applying this approach to nuclear multifragmentation. Then 
the mean number of clans Nc is the mean multiplicity of clusters < M >— < > and the mean number of 
members per clan 

<n> ^ <n>_ 
Nc <M> ^ ' 

is the mean size of clusters. The use of a underlying cluster picture to describe the negative binomial multiplicity 
distribution can be found in Ref. ||39t] . In such a picture the observed particles arise from production of A/ = c clusters 
with probability distribution Pc- This is sequentially followed by each cluster decaying into ka particles with the 
probability Pk^ with a = 1, 2, • • •, c. The probability of observing n = knk = particles is then obtained 

by a compound probability expression 

c 

A negative binomial distribution can be obtained when Pc =< c >^ e^^'^^/c! and Pk^ = {q^" /ka)/hi{l/p). Here 
c = M and < c >=< M >, p = {I - z) ^ (1+ < n > /x)^'^ and q ^ I - p = z = (< n > /x)/{l+ < n > /x). 
Also Nc —< c >— a;ln(l+ < n > /x) and Uc =< n > /Nc ^ {< n > /a;)/ln(l+ < n > /x). This structure can be 
generalized as follow. 

Since the generator of Poisson distribution is an exponential, i.e., the expansion of exponetial gives the Poisson's 
distribution P^{n), 

= E ^ E ^-^^ - E ^r.(AA) (102) 

The grand partition function or the generating function Z = for any distribution can be represented as a Poisson 
distribution whose mean value is the void variable /o or the grand potential 17 = ^/o- the other hand and in 
general we can rewrite the void variable as 

/o(f) = lnZ(x) = ^Xfe=AA^n-(af) (103) 



fc 



where N — ~ X^fc ^fc 



Vk{x) = ^ (104) 



The Vkix) can be connected to its generating function Gix, u): 
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(105) 



Thus the generating function G{x, u) of P„ can be expanded in terms of Vk as 

G{x,u) = ^(1 - urP,,{x) ^ -L^ ^Z„(x)(l - uT 



A/=0 



1 wg] 



M 



M 



M 



E^M(AA)$]n[(i-")'^'»r 



(106) 



M 



Hm k 



where P^.j{J\f) — e ■^A/'^^/M! is the Poisson distribution with the mean oi Af — fa- Here the sum over hm is the sum 
over partitions ft with a fixed M = Uk- Thus we have 



Km k 



(107) 



M 



with M — '^i^Tik and n = fcrife. Any distribution obtained from a generating function of the form of Z{x) = 
e-f" = exp[^j, Xk] can therefore be decomposed as a compound Poisson's distribution with some other distribution 
Vk = Xk/fo obtained from the weight Xk- 

The sequential nature of a process is expUcitly shown on Eg. ( 107 ). The observed particle multiplicity distribution 
arises from a two step process in which M = Uk clusters are first distributed according to a Poisson distribution. 
This is then sequentially followed by breaking each of the nk clusters of type k into k particles with probability 
Vk = XkjN and with n — ^^fc- "^^^ probability associated with a given M and n with x is PM{x,n) — 

PLmUkK" = PLiJ^mki^k/M)-^. 

As an illustration we consider the LC model with Xk — xCkz'^ 12^^^ Using the evolutionary variables |3|] 
X = f3/Ap and z = Ap{l - p) then A/" = a^fc /? for p < 1/2 as already noted so that = e'l^fi^ /M\. The 
Xk = pCkp'^^H}- — p)'^ so that Vk = Xk/N — Ckp''~^{l — p)''- The underlying diagram associated with Vk are 

Xk 



XZ' 



'■ jk and thus Vk = Xk/M is generated from 



shown in Fig.||. For a negative binomial (NB) distribution, 
M = Xk = —X ln(l — z). Therefore the NB is a compound Poisson-Logarithmic distribution as shown in Table III. 

As another example, we consider the HGa model with general a instead of a = 1/2 for LC or a = 1 for NB. The 
weight Xk has the structure of the probability Pk{x, z) of NB distribution given by Eq.(p7|), i.e., 



Xk 



z^V{a + k- 1) 



1 



fc! 



T{a) 



a - 1 (1 - 2) 



— Pr{a-l,z) 



(108) 



thus Vk — [1 — (1 — z)°~^]~^P^^ for HGa. Therefore the HGa P„ distribution is a compound Poisson-NB distribution. 
This may interpreted as a sequential process in which clusters with a Poisson cluster distribution Pc breakup into 
particles with a particle distribution Vk given by a NB distribution. For the various models considered here with their 
Xk listed in Table |, the corresponding distribution Vk and the normalization factor Af — fo are listed in Table III . We 
can further see that HGa can be looked as a compound Poisson-Poisson-Logarithmic distribution, i.e., a distribution 
having three sequential steps; Poissonian breakup into clusters — s- Poissonian breakup of each cluster — > logarithmic 
breakup of each of them. 



TABLE III. Poissonian sequential distribution for various models of Table 



Model 


Weight Af 


P 


<n> 


Geo 




NB 




LC 


2x [1 ^ Vl - z] 


HGa 









Distribution Vh 



JV 



Comments on Pk 



{l-z)z'^-' 



1/2 



z" r(fc-i/2) 



1- VT^ 



r(i/2) 



(l-z)l 



-az z 



Monomer only or Uniform for A'^ species 
Uniform with constituents 
logarithmic with constituents 
NB with constituents with Po = Q 



NB with constituents without k 
Poisson (exponential) 
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Because of a unique role played by the Poisson distribution and the form = exp[^^ Xk\ of the generating 
function, the cluster distribution is usually taken to be a Poisson. However, as noted before, other divisions are 
possible. Using the same approach used above for Pc = Pfj, we can expand any distribution using a NB instead of 



Poisson, i.e., with Pc = P^^ using the form of the generating function Z{x, z) = (1 
by M{z), the normalization factor of a new distribution, we have 

Z{x, z) = [1 - M{z)]-^ = [1 - M{z)g{u = 0)]" 



The g{u) is 



while the G{u) is 



AA([1 - u]z) ^ 



N{z) 



G{u) = j2{^-urp^ 



Z(x,\\~u\z) [l-AAg(u)]~" 



z) 



for the NB. Replacing z 
(109) 

(110) 
(111) 



Thus we have 



M=0 



M\ T{x) 



[g{u)r= Y.pT{xM)[G{u)\ 



M 



M=0 



M=0 



NB 
M 



]=0 



M 



{"felAf fc 



M=0 



p„(f)=^pr(x,AA) ^ n^r(^) 

{nk}M k 



(112) 



(113) 



The result of Eq.(hl3h shows that the distribution P„ can be written as a compound probability distribution of a 



negative binomial Pj^j with another probability Vk distribution generated from g{u). For the case oi J\f{z) = e^, 
which may be considered as the fugacity = e'' for a particle with chemical potential fi ^ z, Q can further be 
decomposed as 



Af{{l-u)z) e(i-")^ 



AA(z) 



= ^(l-«)^pP(z) 



(114) 



fc=0 



i.e., Vk for this case is Poisson P^{z). li M = 1 — e-^"/^, then P^ = Pkifo/^) without fc = and Z{x,z) = 
[l - (1 - e^°l^)Y = ef°. For /o = E J, Xfc given in Table |, the Vk becomes the same probability with Xk replaced by 
Xk/x. As an example the HGa with Z{x^ z) — exp ^]-r^[l ^ (1 ^ ^)^ can be decomposed as a sequential process 

consisting of a NB distribution of clusters with Z{x,N = 1 — eJ°/^) which is then followed by a breakup of clusters 
distributed with a HGa distribution with Z{x = l,z) but without voids, i.e., with Vo = 0. Thus this decomposition 
separates the parameter x assigned to cluster from other parameters. 



D. Poisson transformation and other transformation 



Compound distribution such as those considered in the previous section can also be understood using the fact that 
the Laplace transform of G{u) is related to the Poisson transform for the probability distribution P„ |39|; 



G{un)= dyf{y)e-^y^^ dy f {y)G'' {~uyn) (115) 
Jo Jo 

where G^{n) — e~"e("'^~"^" is the generating function for Poisson's distribution and 
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Jo 



/ dyfiy) = 1 
Jo 

dyyfiy) = i 



(116) 



(n — m)! 



For a Poisson distribution, 



f{y) = 5{y-l) 



G^iun) = £(1 - urP^in) = £ liL^e" 



(117) 
(118) 

(119) 



Tl = 



n=0 



By a Laplace transform or Poisson transform, 



oo />C30 /-oo 

G(un) = V(l - w)"P„(n) = / dyf{y)e--y'' = / dy/(y)C?P(wyn) 

n=0 "^0 "^0 

= ^(l-</ dyf{yy^e-y^ = J2i^-ur dyf{y)P^{yn) 

P„(n)=/ dy/(y)i^e-^«= / d2//(2/)pP(2/n) 
Jo io 



(120) 



(121) 



where P^{yn) is the Poissonian probabiUty with a mean of yn. The Poisson transform of f{y) is Pn{n) and the Laplace 
transform of f{y) is G(un). Thus /(y) is an inverse Poisson transform of Pn{n) or an inverse Laplace transform of 
G{un). However the probability Pn{n) which is a Poisson transform of f{y) may also be considered as a Laplace 
transform of f{y) ^^^^ instead of f{y) itself. The probability P„(n) is a superposition of a Poisson distribution 
P^{yn) with weight of f{y). Thus we may interprete the Pn{n) with a mean number n as the probability in an 
ensemble of mixed systems with various values of the energy or temperature which is distributed with weight f{y)- 
Each of the system with a fixed energy or temperature breakups to give a Poisson distribution P^{yn) with a scaled 
mean number yn. 

Since the G{z,u) = Z{{1 — u)z)/Z{z), we have G{z, w = 1) = Zo{z)/Z{z) = Z~^{z) for the case oi Zq = 1 and thus 



Z{z), 



/ dyf{y)e-y'' = / dyf{y) 
Jo Jo 



Z^{yn) 



(122) 



Thus f{y) is also an inverse Laplace transform of 1/Z{z) with u = l, i.e., 1/Z^{yn) = e In general, the Laplace 
transform and its inversion are 



/•oo 

F{s) = / e- 
Jo 



'f{t)dt 



e''F{s)ds 



Some examples are 



F{s) 



(s+a)" 



2V7rt3 



fit) . 

'("-1)^ 

exp 



constraint and comments 



At 



(n= 1,2,3,---) 
{k > 0) LC 



NB 



(123) 



(124) 



For NB with Z'^"(x, ^) = (1 - z)'"' and n 



G{x, z, u) 



(1 - zf 
(1- (l-u)^)^ 



1 + M- 



z-^{x,z) = {1-zy 



{x + ny 



1-z 
Gix,z,l) 



{un + x)" 



(125) 
(126) 
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Thus is the Poisson transform of 



f{y) = — = X 



[x-iy. 



(x-1)! 



(127) 



which is the inverse Laplace transform of ^"(5 + a) " with s = n and a ^ x and n = x from Eq.( |l24| ) as shown in 
Ref. Q. W ith X = 1, the NB becomes Bose-Einstein (BE) distribution with f{y) = e"^. For LC with Z^'^{x,z) = 
g2x(i-vi-z) and n — — , 

G'(x, z, u) — exp 



2x^J\ - z 1 - ./l 



X\/\ — Z \/\ ~ z 



e^^^^exp 



l\l x\f\^z\l x\f\ — z + 
Z^^{x, z) — exp [— 22;(1 — Vl — z)] = exp 



(128) 



2x^/l - z 1 - 



1 - z + z 



e^-^exp 



1 - z 
G(a;,z,l) 



Thus P'"'^ is the Poisson transform of 



f{y) = 



/T^V^^Vl - z 



exp 



yxVl — z 



(129) 



(130) 



which is the inverse Laplace transform of g-fcv^g2a;Vi-« -vvith s — n + xa/I — z and fc = 2\/ xvT"— z from Eg. (124). 

Similarly we can also decompose any distribution as a superposition of NB distribution instead of superposition 
of Poisson distribution by replacing G^{yn) = e~^" in Laplace transform by the generating function G^^{x,z,u) = 

(1 + "izrj) — (^ x+un ) distribution. That is using McUin transform F{s) — f{y)y^~^dy. 



n=0 
/•oo 



G^^(un,x) = ^(1 - w)"Pr(",2;) = — ^ = X- {un + x) 

71—0 

G{un, x) = ^(1 - w)"P„(n, x) = / /(y')a;"(V")^""+'^"'t^y' 
f{y + x/uh)x^{uyn + x)''^^^^^^^dy 
f{y + X /un)G^^ (uyn, x)dy 
P„(n,x)=/ f{y + x/un)P^^{yn,x)dy 



(-x+l)-l 



(131) 



(132) 
(133) 



with s = 1 — X. Thus a NB transform for the probability P„ may be defined as a shifted Mellin transform for the 
generator G(it). Further investigation of these transformation properties may give interesting features of various 
distributions. As a special case of Mellin transform (s - 1)! = e-H'-^dt. 



V. CONCLUSION 



Event-by-event studies from ultrarelativistic heavy ion collisions such as done at RHIC are being used to study 
the details of particle multiplicity distributions as, for example, those associated with pions. Such studies not only 
give information about the mean number of particles produced, but also information about fluctuations and higher 
order moments of the probability distribution which are important tools for studying the underlying processes and 
mechanisms that operate. They are also useful in distinguishing various phenomenological models. Issues associated 
with fluctuations play an important role in many areas of physics and departures from Poisson statistics are of current 
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interest. One purpose of this paper was an investigation of various models of particle multiplicity distributions that 
can be used in event-by-event analysis. This study was done using a generalized model based on a hypergeometric 
series (HGa) and uses a grand canonical ensemble as its basic framework. This framework has its origin in a Feynman 
path integral approach and involves a cycle class decomposition of the permutation symmetries that originate from the 
underlying wavefunctions associated with the produced particles. Many of the existing distributions used in particle 
phenomenology are shown to be special cases of this more general hypergeometric model. Several quantities are 
shown to be quite general; such as heirarchical scaling relations (Eq.([70|)), initially discused in terms of a particular 
distribution such as the negative binomial distribution. Various models and associated distribution can be developed 
in a unified way. These include the Poisson distribution coming from coherent emission, chaotic emission producing 
a negative binomial distribution, combinations of coherent and chaotic process leading to signal/noise distributions 
and field emission from Lorentzian line shapes producing the Lorentz/Catalan distribution which are all shown to be 
special cases of the HGa model. The HGa model and its associated special cases are used to explore a wide variety 
of phenomena. These include; linked pair approximations leading to heirarchical scaling relations on the reduced 
cumulant level, generalized void scaling relations, clan variable descriptions and their connections with stochastic 
variables associated with branching processes, KNO scaling behavior, enhanced non-Poissonian fluctuations. 

Our results show that even though various distributions have the same mean and fluctuation the distribution itself 
or the underlying mechanism could be very different. Thus to flnd the correct distribution and underlying mechanism 
from the data more informations than just the mean and its fluctuation are necessary. The model used here, HGa, 
has three parameters, x, z, and a. Thus with a given value of a, there is no extra controlling parameter for a given 
mean and fluctuation. All the higher moments and cumulants are then determined by the model without any further 
controling parameter. Further studies of general models having more parameters are needed where the higher moments 
can be used as extra conditions in comparing various distributions. 

We compared various pionic distribution within a generalized hypergeometric model (HGa) which is a special case of 
much more general distribution in grand canonical ensemble. In our model used here for a grand canonical ensemble, 
the weight factor Xk for each species is the same as the mean number < Uk > of the species in the grand canonical 
ensemble. Thus we may determine the weight Xk through experimental data. This result also shows that any power 
law behavior in < > is directly related to the power law of Xk- 

It is explicitly shown that KNO scaling works only for large fluctuations for all the distribution related to the HGa 
model. Comparison within the HGa model also shows that just comparing void variables x ^n.^ ^ is not enough 
to distinguish different models that describing pion data. Thus new parameters have to be found which are quite 
different between different models. Beside mean values and fluctuations higher order reduced factorial cumulants need 
also to be evaluated. 

In this paper we have also generalized the Poisson transformation and the compound distribution that arises 
from sequential process. Specifically, the underlying sequential picture involves a two step process where the final 
distribution arises from the production of clusters followed by a subsequent decay of the clusters. For the HGa model, 
the final distribution is obtained from compounding a Poisson distribution of clusters with a NB distribution coming 
from the decay of each of the clusters. The HGa may arise through a three step sequential process of Poisson-Poisson- 
Logarithmic compound distribution. It is also shown that the HGa can arise from a two step sequential process of a 
NB distribution followed by a new HGa with a different mean value. 

This work was supported in part by Grant No. 2001-1-11100-005-3 from the Basic Research Program of the Korea 
Science and Engineering Foundation and in part by the DOE Grant No. DE-FG02-96ER-40987. 
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